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ABSTRACT 


During the design stage of automatic control systems for 
aircraft, fire control or other physical systems, system com- 
ponents must be selected. Values of the parameters of these 
components must be determined so as to produce the desired 


dynamic properties of the resultant system. 


The method of estimating dynamic characteristics employ- 
ed in this investigation is based on the encirclement criteria of 
Cauchy's Residue Theorem... A computer was developed to aid in 
this estimation by making certain modifications to the data pre- 


sentation of the Fourier Harmonic Synthesizer. 


This investigation indicated that dynamic characteristics 
may be readily estimated by the method to be described for any 
selected set of system parameters. System parameter specifi- 
cations and limits may also be set when the dynamic character- 


istics are specified. 


Several examples are included to illustrate the techniques 
involved and to indicate the suitability of applying this method to 


aid in the analysis of higher order systems. 
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OBJECT 


The object of this investigation is to presént suitable 
techniques, based on the Cauchy Residue Theorem and utilizing 
a modified Fourier Harmonic Synthesizer as a computer, to aid 


in the estimation of the dynamic characteristics of physical 
systems. 


ll 





CHAPTER 1 


INTRODUCTION 


In the design of physical systems some of the most important 
considerations are the dynamic properties of the resultant system. 
It will be shown in Summary 1 that the differential equation des- 
cribing the system may be written with the input terms on the 
right side. By setting the right side equal to zero the homogeneous 
equation is obtained. This equation determines thetransient modes 
of the system, which are independent of the input. The actual in- 
put affects the magnitude of the transient modes but has no effect 
on the damping ratios, the oscillation frequencies, or response 
times of the transients. These characteristics are determined 
by the component parameters. Once the system is built, control 
over these dynamic properties is restricted. Since these proper- 
ties have a prominent effect on the performance of the system, the 
engineer must have some means of predicting them during the 


design stage. 


Exact determination of system parameters is often neither 
necessary nor desirable in the early stages of design if limiting 
values of the parameters can be estimated interms of the dynamic 
characteristics of the system. Such an estimate will enable the 
engineer to set limits on the component specifications in the in- 


itial stages of design. 


The method of estimation of system dynamic characteristics 
described in this thesis is based on a concept first advanced by 
Sidney Lees in Massachusetts Institute of Technology Instrumen- 
tation Laboratory Report R-71 of April 1954 (1) 


The Mapping Theorem, described in Instrument Engineer- 


ing Volume II ©) and in Theory of ServoMechanisms, (3) and 
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the Cauchy Residue Theorem described by Guillemin (4) and 


others, are the basis for this concept. 


Lees applies this theorem (1) and describes a procedure by 
which dynamic characteristics may be estimated from polynomial 
contours approximated by the computation of certain selected 


points. 


If these contours can be accurately computed and displayed 
by some device, the application of the Cauchy Theorem to this 


problem should be greatly facilitated. 


A Fourier Harmonic Synthesizer built at the Massachusetts 
Institute of Technology Instrumentation Laboratory and described 
by Instrumentation Laboratory Engineering Memo 6445-E-42 (9) 
was used for this purpose. Modifications to the existing instal- 
lation and the techniques employed to utilize this device as an 
aid in the estimation of the dynamic characteristics of physical 


systems are described in the succeeding chapters. 
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CHAPTER 2 


MATHEMATICAL BACKGROUND 


The performance of physical systems may often be described 
by linear differential equations with coefficients that are constant 
once they have been decided upon. The properties of the compo- 


nents of the system may be obtained from the equation coefficients. 


The complete solution of the linear differential equation 
consists of the transient solution, plus the forced solution as is 


shown in Summary l. 


The transient solution of the linear differential equation is 
the solution of the homogeneous equation, which is obtained by 


setting the right hand side of the differential equation equal to zero. 


SUMMARY l 


As expressed in Instrument Engineering Volume II 5) , a 


typical ordinary linear differential equation has the form 





n i m k 

a =~ = », b, <E (1) 
ES, dt dt 

i1=0 k=0 


The homogeneous equation is 





n 
i 
) a, =~ = 0 (2) 
a dt 
i=0 


The associated characteristic equation is 


ar z+ = 0 . (3) 
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A fundamental theorem of algebra described in Summary 
24-1 Instrument Engineering Vol. II (2) states that every 
equation of nth degree has n and only n roots. The roots of 
Equation (3) are Zi, Zo, Zags +++. By. In general the roots are 
complex numbers so that 
Z =X, + Ua (4) 


1 


When all roots are distinct, transient solutions of Equation 


(1) have the form 
is Zt i‘ i : 
Veer) 2) cre * = » ce 1 e (5) 
“s Ls 


1=1 i=l 


Roots of the characteristic equation must have the di- 
mensions of inverse time (from Equation (5) ). Real roots may 
be expressed as 


2 ae 
22%, 2 -—n = <2 (5) 
1 


where (CT); ae characteristic time of real root. 


Complex roots may have the form 


| F 2 
i =e Wi = = Se Te) oe (7) 


kk a's 
where OK (DR), = damping ratio of kth root, 


os undamped natural frequency of kth root, 
k 
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- a eee S = 
and x = - o) in = - where (CT), ie 


characteristic time of kth root. 


In polar coordinates, the kth complex root form is 


ee ™- (ADR), | 
"Xk 


) 
ok 8 
ca g (8) 


where @ = Ff = (ADR), = 7 = o 
and (ADR), = A, = cos’' {, = angle of damping ratio of the 
k 


kth root. 


The transient solution determines what oscillatory modes 
can be present, and their mode shapes, and also determines the 
extent to which the dynamic performance of the system will be 


affected by the specified input. 


The forced solution reflects the effect of the input and 
determines the relative magnitude of the modes indicated in the 
transient solution superimposed on the forced solution. It also 
determines the system response time due to the effect of the in- 
put on the mode shape but has no effect on the damping or oscil- 
lation frequencies of the second order modes. The response time 
for the system will, in no case, be greater than that indicated by 


the slowest decaying term of the transient solution. 


This investigation is not restricted to any specified input 
and considers only dynamic effects resulting from the roots of the 
characteristic equation. The analysis of this part of the complete 
performance equation gives information on dynamic properties 


which greatly affect the performance of the system for any 
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arbitrary input and it is to this end that the investigation is 


directed. 


The roots of the homogeneous equation expressed as a 
characteristic equation with the variable z are defined in M.1.T. 
Instrumentation Laboratory Report R-7l (1) as the dynamic 
characteristics of the associated system. The dynamic charac- 
teristics can be determined if the roots of the characteristic 


equation can be found by factoring or other means. 


Lees (discusses the important considerations in the 
investigation of transient solutions of most practical systems. 
One consideration is that the transient solutions decay to zero, 
thus providing absolute stability. In addition, it is necessary 
that the transient solutions decay rapidly enough to meet specified 
conditions; this is the relative transient stability requirement. 
The damping ratio and the oscillation frequency are also impor- 
tant considerations in investigating the transient solutions. 
These are the dynamic characteristics and are apparent from the 
transient solution or roots of the characteristic equation. In this 
thesis, the terms 7, © and We specified by the roots of the charac- 
teristic equation are considered to be the dynamic characteristics 
of the physical system. The equivalence of these terms and the 


characteristic equation roots are shown in Figure 1. 


Imaginary 
Axis (y) 
Real Axis (x) 





Real Root Complex Roots 


Figure 1; Roots of the Characteristic Equation as Represented 


on the Complex Plane, 


Absolute transient stability of the system exists when all 


of the roots of the characteristic equation have negative real parts, 
: xt x,t 
i.e., when x, and xX, are negative, e and e décayto.zeroas 


t — 0. 


Relative transient stability exists when the smallest of the 
magnitudes of the real parts of the roots is greater thana 


! 
specified reference magnitude, i.e., if x, * | > | then 


mp 











x 
ref 


(CT); < (CT) ce (%; is most positive real part of all the roots). 
mp 


This has been discussed in greater detail in Instrument Engineering, 
Volume J] Jey) Chapter 20. 


The damping ratio and oscillation frequency of the system 
are related to the real and imaginary values of the complex roots 


as shown in Equation (7). 


In many practical systems the dynamic characteristics are 
restricted to a certain desired sector or range depending on the 
relative transient stability requirements and specified undamped 
natural frequency and damping ratio limitations. These sectors 
and ranges are shown in Figure 2 below. Specification of a sector 
imposes limitations on the roots and thus the coefficients of the 
characteristic equation, and it is desired to determine the range 
of each coefficient which will provide dynamic characteristics 
within specified limits. Figure 2 shows how these limitations may 
be investigated with respect to the characteristic equation complex 
plane root loci. These contours are useful when investigating for 
relative stability, oscillation frequency and damping ratio. In 
most practical systems the specifications impose limitations on 
two or all three of these characteristics and a contour similar to 


Figure 2-d is used. 
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Imaginary 
Axis (y) 


CT ref. 













(a) (b) (c) 
Relative Undamped Damping SuperimpoSing Figures 
Transient Natural Ratio a. b. andc 
Stability Frequency Limits Dynamic Characteristic 
Limit Limits Limits 


Figure 2: Standard Contours 


This investigation is concerned with the techniques con- 
cerning the solution of two problems: determination or estimation 
of the dynamic characteristics, and fixing the limitations and re- 
lationships of the coefficients. Fixation of the coefficients of the 
characteristic equation determines the limits of the component 


parameters in the actual system. 


The simpler case, that of locating the roots of the charac~ 
teristic equation with known constant coefficients will be taken up 
first. This case will be treated in some detail since all the 
principles involved are equally applicable in estimating the dy- 


namic characteristics from equations with adjustable coefficients. 


In this presentation, the methods used in estimating the dy- 
namic characteristics of physical systems are based on Cauchy's 


Residue Theorem, as presented by Guillemin (4) and others, 


From Equation (3), 
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This is expanded as 


n n-1 2 
P(z)=a,zt+a,_)z +...ta,9Z +a, z2+a, = 0 (2) 


From Cauchy's Residue Theorem, if P(z) is a polynomial 
of z, P' (z) is the derivative of the polynomial with respect to 
z, and C is a closed countour on the z complex plane containing 
N roots of P(z) =0, 


a | 
FO) az =520N (10) 
P (z) 
Cc ’ 
This integral may be interpreted to mean that the argument 
of P(z) increases exactly 27N times as the closed contour C is 
traced out in a positive sense, orthat.the P(z}plane path encircles 
its origin N times. This is also referred toas the mapping theorem, 


in Instrument Engineering, Volume II (2) and is illustrated in 







Figure 3. 
Closed contour P(z}plane path 
C in z-plane corresponding 
(a) Three roots on z-plane (b) Three encirclements of the 
within contour origin on the P(z)-plane 


Figure 3: Illustration of Mapping Theorem 
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A closed contour similar to Figure 2-d and corresponding 
to Figure 3~a may be drawn on the characteristic equation root 
complex plane, in which it is desired to investigate for roots of 
the characteristic equation. By substituting values of the 
variable lying along the complex plane contour of Figure 2-d, a 
corresponding contour on the P (z) plane corresponding to Figure 
3-b may be determined. The number of encirclements of the 
P(z) plane origin is the number of polynomial roots within the 
root ‘complex plane contour. As will be shown in the examples, 
the polynomial P(z) of Equation (9) usually includes a constant 
term ay: It is convenient to write the polynomial as P(z)-a, 


and 5 becomes the P(z) plane encirclement point. 


The P(z) plane path of Figure 3-~b may be calculated by 
hand by substituting values of the variable along the contour of 
Figure 2-d. This becomes impractical since a large number of 
calculations is required if a complete and accurate P(z) plane 
path corresponding to the z plane contour of Figure 2-d is desired. 
Chapter 3 describes how the Fourier Synthesizer may be used to 


obtain the information provided by the P(z) plane path. 
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CHAPTER 3 


APPLICATIONS OF THE FOURIER HARMONIC SYNTHESIZER 


The Fourier Harmonic Synthesizer may be used to facili- 
tate the application of Cauchy's Residue Theorem to this problem, 
reduce the required calculations, increase the accuracy and bring 


other related items of interest to light. 


It has been shown that the polynomial in z may be written 
in the form 


(7) = By 3 re | ee Ro ef 100 


jg (11) 


by ‘setting z = Re jo 

This polar form of the equation is suitable for setting into 
the harmonic synthesizer. The method of scaling to obtain full 
scale deflection of the synthesizer and maximum accuracy there- 


by is given and explained in Example l. 


Each of the coefficients (R” ; Ro) corresponds to the 
radius of the contribution of the particular term of each order as 


the variation is traced out from @ = Oto 90=27,. 


The synthesizer traces the polynomial | P (z)|plane from 
6=0 to 96=27 for the selected R. The plot as recorded on the 
Sanborn Recorder is displayed as separate imaginary (y) and real 
(x) components. The polynomial contour (P(z) plane) which may 
be constructed from the . and y components, is symmetrical 
about the real axis. The real component is a mirror image 
about 180° and the imaginary component crosses through zero 


at 180° with the 180° to 360° portion being an iiverted mirror 
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image. Therefore, as is indicated in Figure 4, only the 0° to 


180° BPCrUlon is requir car 





0 cs 27 
Imag. 
Y 0 
Real 
0 
aes 
0 q 27 
(a) Synthesizer Component Traces (b) Polynomial Contour 
Represented by 
Figure 4 Synthesizer Traces 


The P(z) plane may be drawn from the imaginary and real 
components, and the number of encirclements of the constant 
term counted. This, as Cauchy's Residue Theorem states, gives 
the number of roots that may be found at smaller radii (R) than 
that selected. 


The major item of interest is the y (imaginary) cross-over 
(where the y component equals zero) in relation to the constant 
term a, on the negative real axis. This point is represented by 
a horizontal line at x= - a, on the x (real) component plot. Only 
the @ at y cross-overs need be noted from the imaginary com- 
ponent plot, and from the real component plot only the fact of 
whether the corresponding points are more negative or less ne- 
gative than the constant term need be noted in order to sketch 
rapidly and adequately the approximate P(z) plane. If ata y 
(imaginary) cross-over, the x (real) component crosses the con~ 
stant term line at that point, roots occur at that radius, one real 
root if @=7, two complex conjugate roots if 8 equals anything other 
than 7. This immediately gives the (ADR) and (DR) of those roots, 
since (ADR) = 7-@ and (DR) = cos’ (7 - @), as well as the mag- 
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nitude of the roots, since R=w, . The typical P (z) planes of 


Figure 5 will illustrate this. 
, 





0 
sare) 
Three encirclements of a, One encirclement of ay 
3 roots atwi <i lL root at w <8 


2 complex roots at Wt: 8 


6 = 3, C= cos 4 =. 707 
(a) P(z) contour for R=1 (b) P(z) contour for R= .8 
Figure 5: Typical P(z) Contours 
The approximate P(z) plane need not be sketched at all 


since the encirclement criteria is immediately apparent from the 
x and y component traces directly. If during the interval between 
two y cross ~ overs, the real plot crosses the constant line, and 
remains on the opposite side when the y plot again reaches zero, 
One encirclement is indicated. That is, one root on the complex 
plane exists at a radius less than that for which the plot was made. 


This is illustrated in Figure 6. 





0 1 
‘ 
0 
x 
0 
“a 
O 
0 oy G5 T 


One encirclement between O and Oo) 


No encirclements between 0, and 


G6 Ore 
(a) One Encirclement Indicated by (b) One Encirclement In - 
Synthesizer Traces dicated by P(z) Contour 


Figure 6 
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It should be noted that from the component or P (z) plot of 
just a single radius, no conclusion should be made as to whether 
the roots indicated are real or complex, nor estimate the (DR) 
of such roots unless some additional information is known about 
the equation being tested or unless the x plot crosses the constant 
line at or very near the @ at which the y cross-over occurs. 
Caution should be exercised when attempting to draw concrete 
conclusions in estimating (DR) from plots of x plots crossing the 
constant line near a y cross-over until the effects of R variation 
for the particular equation is fully appreciated and reiative mag- 
nitude of synthesizer indications and tape scaling factors am consi- 
dered. Examples (1) and (2) demonstrate, among other things, 
how an error in estimating (DR) may occur by considering only 


the results of a single run. 


By making several runs of constant R for varying 6 on the 
machine, a composite plot may be constructed as is shown in 


Bxample (3). 


The magnitude of the roots may be estimated within the 
limits of the selected R's by the criteria described above. Curves 
connecting points of the various R plots for equal angles (@) pro-~-. 
duce curves that are actually constant 0 curves for varying R. 

The number of encirclements of the constant term by these con - 
stant 9 curves give the number of roots at 6's greater than the 
value of 6 for that curve. This gives the limits on the damping 


ratio, since € =cos (180° - 0). 


As is indicated in Figure 7, the number of encirclements of 
the constant 6 curves is obtained by considering @ to be held constant 
until an infinite radius (R) is reached then permitting @ to complete 
its sweep to 180°. The latter part of the sweep ( @ infinite R) will 
never appear on any trace or plot but must be considered in order 


to determine encirlemenis. 
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-  9=135° 0<R <0 
135° < @<360° at R=a 


(a) Constant 6 curves for 
Third Order Equation 





(b) Constant @ curves for 


Fourth Order Equation 


Figure 7: Typical Constant @- Variable R curves 
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CHAPTER 4 


TECHNIQUES INVOLVED IN THE TREATMENT OF SPECIFIC 
EXAMPLES 


EXAMPLE 1: Constant Coefficient Case 
In order to scale the problem in a manner suitable for. 
setting into the synthesizer, the polynomial in z is written in 


the polar form as shown below. 


The polar form coefficients are determined for each of the 
trial radii (R's) selected. The scale setting is then determined 
so that the largest coefficient may correspond to the largest 
synthesizer scale.setting (1. 0) in order that the maximum syn- 
thesizer accuracy may be obtained. This simple set df slide 
rule ratios lends itself to tabular presentation as is indicated 
in Table I, page 32 which follows the examples. 


P (z) - 52=2° +13 2° + 60 z° + 120 2° + 124 z 
Ree ts Ree > eo: elegy 
+120 RR? e J29 4 124 Re J 


Since z= Re 
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EXAMPLE 1 


Constant Coefficient Case 


The Sanborn record of the synthesizer component traces are sketched and summarized below for convenience. 
The values at the important points and cross-overs are included and the corresponding P(z) plots sketched in order 
that the effect of R variation may be readily observed. The conclusions and estimates based on the individual runs 


are given with the corresponding traces in order that the technique used in the investigation may be followed con- 


veniently. 


P(z) - 52 = 2° +13 2° + 60 z° + 120 2° + 1242 


0° 47° 92° 127° 148° 


5 encirclements; hence 5 roots at R< 6 
y 
At least 1 real root. O 
One complex pair may be expected with 
a 6 at F R=6 
148° > @ > 127 
or x 
127 > 8 > 02° O 
“RY 


2 roots @ R = V26 





o OF Sie lOO®° «145° or 


6 = 145 A, = 35 
; y 
(The run for R=6 indicated the possibility 0 
of roots @ 6 = 145° The other indication, i.e. 
127 @> 92°did not materialize) RADE 
3 encirclements 
3 roots @ R< V26 © 
Due to the position of the loop in the P(z) -52 


plot the only safe estimate is that if there 





is another complex pair, it should occur 
at 100k e<z 


Obs Shale We 7 
i 3 roots @R<3.5 y 
| If there is another complex pair at R< 3.5 0 
it is more likely to occur between 11° and x R=35 
Boe then between 59° and 111°. (The basis of this 3 
me 0° Xx estimate is the position of the constant term x 
and observing the behavior for many other 0 
-52 equations) “oc 
0° S32 50~ 7 


y 

y 

2 roots @ R= V2 0 
= o 2 9 

= 6 = 130 A 2 50 RAD 

130° 0° - l root (real) @ R < V2 : 
e) 
_- 52 
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Example 1 (cont. ) 








0° 118° v 
y y 
0 
18 x 1 root @ R<1.2 Rel.2 
O° 
-62 (The polynomial contour plot and 
component tracea indicate the % 
real root ahould be very close to =52 
R =1, 2) 
o° v 
y 
0 
R=/0 
l root @ R=1 
@=x A,=0 ‘ 
- O 
All roota have been located “S52 
o° v 
y y 
0 
R25 
x ° 
0° No roots occur 
-52 at R<.5 ‘ 
4 0 
Thia checks with previous runa 52 
since all roota have been located 
at R'a > .5 
oe 4 3 2 
The equation P(z)=2z* + 13z + 60z° + 120z° + 124z + 52 
has roots of z=-Stj zz- l+j z=-l 
(wo, = V26, A, 1°) (w= N2,A, = 45°) (@=1,A 70) 


The actual values of the roots of the equation agree very well with those 
indicated by the aynthesizer. 
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EXAMPLE 2 


Constant Coefficient Case 


EE 


Example 2 is of a form very similar to Example l, even to the extent of an increasing magnitude of the 


coefficients with decreasing powers of the variable. 


Example 1; P (z) = 


Example 2: y : 


Boas? + 60 2° + 120 o> + lod eee 


P(z) sz +52' +362" +88 2" + 108 2452 


(The constant term (a, = 52 in both examples ) while greatly affecting the roots has no effect on 
the shape or magnitude of the component traces or polynomial contour plot) 


These two examples, while similar in form, show quite different behavior and indicate that an estimate 


of (or A "4 ) based on a single run may be considerably in error. 


P(z) - 52 = 2° + $24 + 362° + 88 2° 


5 encirclements,5 roots @ R<7 
If complex roots exist, one pair 
may be expected @ 142> > uo 
11d?>@ >84° 


or 


2 roots @ R= V 26 

g=101° A = 79° 
Second estimate from the R=7 run 
is aubstantiated. (Compare the 
R= V26 runs of Example | and 
Example 2) 





3 roots @ R < V26 


3 roots @ R<3 


If more complex roots exist, it 
cannot be stated positively whether 
they will be at 672@<136° 

136% <7 

The relative magnitude of the trace 
peaka indicate however, that 6 will 
be closer to 136° than either 67° or 7 


or 


2 roots @ R= V2. 


@= 135° A, = 45° 
. 4 


1 root (res!) @ R < V2 
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+ 108 z 


45° 84° 110° 142°) ox 


0° 52°, 101° 140°: on 
y 
O 
RN26 
Kx | 
O 
-52 
o° 67° 136° 


0° 105° 135° 
y 
e) 
R=12 
x 
0 
coe 


Example 2 (cont. ) 


O° 
y 
y 
l root @ R=l O 
6=7 A, = 0 Re| 
0° X% All 5 roots have been located. 
-55 No indication of roots should be =A 
found at R's <1 -52 
Ot 
y 
y 
0 
R=5 
oP x No roots @ R<.5 - 
=52 O 
-52 


The polonomial: 


P(z) = 2° + 5z* + 362° + 88z- +108z +52 


has roots: z=- 1+ 5j z~=z-ltj zz-l 
(e) °o 
(w= N26 , A,=79") (wy V2, A, =45°) ( @ = 1,A,/=0°) 


In this example the actual values of the roots agree almost exactly with 


the values indicated by the synthesizer. 
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EXAMPLE 3 


Constant Coefficient Case 


P (z) - 120 = 2°02? + 442 


y OF 85° 140° 
y 
0 
85 : 3 encirclements R= 
° ° 
0 0 3 roots @ R<7 : 
-120 (Possibility of one real root is 
indicated at R less than but -l20 
close to 7) 
( Suspect one complex pair @ 
85 °< @< 140°) 
y OF S57 loo ee 
y 
0 
95° : R=6 
isa. OF 
1 root (real) @ R=6 x 
2 tao 
120 O=" A, = 0 0 
-120 
2 roots @ R < 6 
If a complex pair exists it should be 
found in the vicinity of @ = 139° : 
y 0 100° 142° aw 
y 
0 
2 t < 5.2 
iQ0° Oot Oe ae R=425 
x 
142° 0° 
x 
=120 
0 
=12O 
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Example 3 (cont. ) 


0° 16° 134° oor 

y y 
0 

NG R=4.47 

x 

134° 0° x 

“120 | os 
2 roots @ R = 4.47 
@= 1169 A= 64° 
5 
O° T 
y 

y 

0 
R= 

0° x x 
-120 Re 
No roots @ 120 

R<3.5 
The polynomial 
P(z) =z> + 1027 + 44z + 120 
has roots z=- 6 z=-2+ 4j 
: fo) 
a ee 0°) (w, = 4.47, A, = 60°) 


The factored form is therefore 
P(s) = (z + 6)(Z + 2 + 4j)(z + 2 - 43) 


The accuracy of@ indicated by the synthesizer is very good except when 6 approaches 180°, 


Figure 8 shows all the polynomial contours sketched above. It was constructed by taking 
the values from the tapes to provide a series of superimposed polynomial contours for all 
values of R, 

Another series of curves was constructed by connecting points of equal @ in order that 
constant @ curves for varying R may be obtained. Cauch's Residue Theorem is applied to 
the constant 6 curves. The encirclements are counted and the @ values determine the sector 
in which roots, if any, are located. 
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EXAMPLE 3 


9=90° Gay 
50 
= ee seal 
roots ; 
a g@>90° “a 40 
ye 9°135" 304 92120° 







| encirclement 
1 root (reol) 


§=m(reol root) 
tim (P(z)-120] 











30) -20s 10 





120-0 100-90 -€ 


Gg -104 Re[P(z)- 120] —= 


r=5.25 50 


P(z)-120 = 23 +1022 +442 = 30/29 +}Orte i244 
44re/8 


100 
sie) 
=|20 
“130 
“140 


“150 


Figure 9: Expanded Portion of Region Near Origin in Figure 8 for Constant @- Variable R. 
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EXAMPLE 4: Adjustable Coefficient Case 


Example 4 illustrates one analytical method of investigating 
a typical equation with adjustable coefficients in order that the 
limits of the component parameters affecting the coefficient may 
be determined which meet the specified dynamic conditions. This 
method as explained in Example 4 utilizes the algebraic relations 
between the adjustable coefficients and the sums and products of 
the resultant roots. This method is suitable for lower-order 


equations (second or third) with one or two adjustable coefficients. 
EXAMPLE 4 


A typical characteristic equation,of a positional servo- 
mechanism with modified velocity signal damping, is: 


27 r 'p 
zea |P haa ——~————— + 2 
P (ps) 1+ 27 r 'p 





t - 
© (ps) (res)} P+! = © oo 
where 'p = = P > PR = Laplace Transform operator 
n 
(ps) - 

So W = undamped natural frequency 

#7 HS) of positional servomechansim 
e (ps) = damping ratio of positional 
servomechanism 
l Ty: 

ar = CT = ae = non-dimensional character- 

eee (ps) a5 (ps) istic time of signal modifier 
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Expanding Equation (1): 


3 1 _@ 
: *|P© (ee) * aca) Tee a A 











L ee 
2 £ (ps) (res) Dar ry 2 ree (2) 
V V 
Or 
1 3 1 2 ! = 
DS laze. b) *p + (be+ 1) pitbp=-= 0 (3) 
where 
=2 
o * *(ps) 
: l 
£ 2 r 
Vv 
c= 2 
> (ps) (res) 
For zero forced dynamic error, 
Za 
ED Ey ene) A - 0 for finite w 
WO (cm) n (ps) 
( ps ) 
where A (cm) (angular velocity of controlled member) + 0 
Thus, ¢ (ps) (res) c= 0 
and Equation (3) becomes 
3 hts i 
Dp 4 Ca + DP) ape ee et = 0) (4) 
Substituting 'p aes EF , letting p = We e J z 
n (v) 
(ps) 


40 


and defining R = — he eae a pay 
n 
(ps) 
; ; 2 3 y) : 
Equation (4) is: P(z)= z° +(at+b) z*~“+z+b=0 (5) 


If the roots of the polynomial in z [P (z)] are 21, Zo, and Za; 


P (z) (2-2)) (2-25) (z-2Z.) = 0 (6) 


3 2 
Z -(Z2,+2Z5t Za) Z + (2, Zo + Zo Zag + 2) 2, )2- 


Z) Zo Zag = 0 (7) 


For a third order equation with two complex roots, the roots 
will have the form: 


By SS ia 
. ; 2 
Bigs GE MIG SS Wn see WS 
ey (vr) 
2 ae , AF 2 
3 7 QO Ue 8 os Oe - JM, tae eve 
(r) (r) 


Using these relations and equating coefficients of like terms in 
Equations (5) and (7) we have: 


amiaUEeent 2) Comine (8) 
(r) 
2 
|, = fe + W (9) 
ren py)”  {r) 
2 
b = (10) 
O “Tiny 
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For the purpose of illustration, the specifications for the 


dynamic characteristics were taken to be 


l) a Sw) <I 


iV 
no 


setting C .. =o ae 


From Equation (8) a+t+b = a@+t+w 








Using Equation (9) l - we aa (atb) - Wo 
(rv) (vr) (r) 
or W = : 
Der) a+b 
Thus @ = atb- : _ (a¢b)*- 
atb . a+ 
Pee (alt bien 1 (a+b) * -1 
" (r) athb (a+b) 2 (a+b) 3 
Choose a+b, calculate b. 
meus ele bind l ae ee (acaba 
niex : (atb)? {a tb) (a+b) 
 (atb)l , - 3 
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and b = = = 
ax (43) ° 3 
Setting w | = ,4 
(r) 
When ree SOs atb=at, 4 
a .=(a+tb) =.4 
1 (4)? = 2 (.5) (4) [ (arb) -. 4] 
atb= S24 04 = 2.5 
@=2.5~.4 = 2.1 
: 2 
b = (2.1) (.4) = .336 
When ¢ = 1. atb= a+.8 


a=(atb)-.8 


lan aes 2 (1) (.4) | (at b)- . 8 


atb = 8 


3 1.85 


@ = 1,85-.8 = 1,05 
2 
b = (1.05) (.4)* = .168 


Other points between these limits are calculated similarly by 


using other values of .5< Gas l 
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For @ = .2 b=aw, = .2u, » Wy = 
(vr) (vr) (vr) 
iol.|6U= 2 ee (. 2) 
n r oh 
(rv) (vr) 
L-0? 
Zag & & W = eee). = AUCs ) 
Cr en ee n 
(r) (vr) 
, 2 
apb=.2 + 20 eu = ,.24 5 (l-w ) 
n n 
(vr) (rv) 
= .2+5(1-5b) 
=72 +5 = 25) 
a+b = 5.2 - 25b 
b 9.2-(artb) 
ao 
Choose values of a+b, calculate b. 
For C . = | atb =at 20, » le we = 2W (atb)-2w , 
eS (rv) (rv) (r) (r) 
_ 2 
b=@ WwW = |(atb) -2 WO, WO 
(rv) r) (r) 
hy 
z (v) 2 
2, Nir) 
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(a W 
Z (vr) 
b W 

: "(r) 

l= w? 1+ 3 we 

a+b = se + 2w Sk ne | 
Yn "(r) 20), 

(vr) (vr) 


Assume 0sw, <1, calculate corresponding a+b and 
1B 


bong. =1 line. 


From the preceding calculations, Figure 10 was constructed, 
which shows the possible combinations and variations of the 
parameters (a+b) and (b) which satisfy the performance specifi - 


cations. 
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The complex plane root loci in Figure ll may be constructed, 


with the variables (atb) and (b), as shown below: 
P (z) = 2° + (a+b) o* +2 +b = 0 


Relating equations: 


ad D> @ + 20 .w, (a) 
(r) 
Ll =wi +200 a (b) 
(r) (r) 
b =a as (c) 
n 


Solving (b) for a: 


2 
l-W 
n 
Go 
gamle 0: 
(r) 
Substituting in (a): 
l-we 
Asya to oe tw (2) 
2 Ce W ron (rv) 
(r) 
Solving for 2¢, WO, 
(r) 
2 2 
(atb) t NV (atb) - 4 (1-0) ) 
r 
32. ee a ee Up Fee (3) 
(r) : 
By assuing 0 < WO <0 ee Ze WO can be calculated 
r 


(r) 


AT 


in terms of atb. 


For each assumed value of w , values of a+b may be 
(r) 
tabulated, and the corresponding a's and b's found which satisfy 


the above relating equations. 


These tables are presented in Appendix B and the results 


are plotted in Figure 11. 


The qa = .2 line was obtained by letting a = .2 in Equation (a). 


Solving for ¢ oe 


(r) 


The points on the q@ = .2 line will lie at the intersection of 


the a+b = constant curves and the corresponding vertical Ce WO 


(vr) 
lines. This curve is also plotted in Figure 11. 


The area enclosed by the a = .2 line, the ¢ = S .9 line, the 


GG. =1.0line, and 4 Lu, is the area in which the values 


(r) 


of a+b and b will satisfy the specifications of the problem. 


r 
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One might also construct such a plot as Figure ll solely 
through the use of the synthesizer. Runs are made at various 
radii, (R), and for various values of the second coefficient, (a+b), 
and reading the values of the constant term, (b), which would give 


roots at each combination of (R) and (atb). 


This is not recommended, however, since a very large num- 
ber of runs would be required, with the number of runs required 
being further increased if there were more adjustable coefficients. 
Analytical methods such as this, or use of the cubic charts in 
Volume II, ) Instrument Engineering would enable one to construct 


a coefficient boundary plot in less time. 


In engineering applications and control system design prob- 
lems such stringent boundaries may not be required. If the pri- 
mary requirement is selection of components with parameters (repre- 
sented by the coefficients) which will permit the system perform- 
ance to be well within the specifications, use of the synthesizer 
may greatly expedite selection of coefficients which appear to be 
satisfactory. Checking this selection on the synthesizer will insure 


that the performance is within specifications. 


In this instance only two runs On the synthesizer, one ata 
radius, (R), of about 1.0 and one at a radius of about 0.7 for an 
(a+b) of around 1.5 will be sufficient to give a rough but adequate 
locus of complex roots on the root plane. Two more runs at an 
(a+b) of about 1. 7 at the same radii will give another locus. These 
four runs and two resultant loci will be sufficient to give a general 
idea of the loci of other selected values of the coefficient. Figure 
12 shows such an estimated plot which corresponds to the computed 


plot of Figure ll. 


90 








Approximate Root Plane Plot of Example 4 


Figure 12 
From this rough root plane locus it is apparent that a value 


of (atb) of 1.6 with a constant term (b), of 0.2 should be one of 
many possible ‘combinations which give a satisfactory solution 
well within specifications. There remains only to set up this 
equation to check the proposed selections. Few trial and error 
runs will be required to locate the roots accurately, since a 
glance at the rough root plane plot will show that the complex pair 
should be found at an Ws» of about 0.65 or 0.7. A rudimentary 
knowledge of Root Locus will indicate that the real root should be 
found at a radius of about 0.5 since for this equation the real root 
moves left from the origin at twice the rate the x component of 
the conjugate pair moves right. The Root Locus method is de- 


scribed in Volume II, Instrument Engineering (2). 


As may be seen from Figure 10 giving the boundary values 
of (a+b) and (b) which yield acceptable solutions, the three adjust- 
able parameters, a,b, andc, reduce to two adjustable coefficients. 
It is shown that boundary values of two adjustable coefficients may 
be presented on a two-axis or two-dimensional plot. It follows that 
with three adjustable coefficients, which may be made up of any 
number of adjustable system parameters, a cubic or three-dimen- 
sional plot would be required. 


ol 


It is apparent that with more than two or three adjustable 
coefficients, the presentation of definitive boundary conditions 
for each coefficient when dependent upon each of the others be- 
comes impractical. It was also found that a method based on 
the algebraic relation between adjustable coefficients and re- 
sultant roots becomes unwieldy for fourth order or higher 


equations. 


I 
eo 


i.e., 158t coefficient 


9 nd coefficient 


tl 
ISA “3 

N 

p- 


n 
nth coefficient = us Z 


Accordingly, it appears that the most generally suitable 
method of investigating the result of adjustable parameters, and 
hence adjustable coefficients, will be to set up a series of in - 
equalities for various conditions and observing the effect upon 
the resultant performance. From these inequalities reasonable 
engineering estimates may be made for the selection of various 
sets of coefficients, and the resultant dynamic characteristics 
that each set should yield. These estimates are readily checked 
on the machine in the manner given for the fixed coefficient case 
to insure that they do in fact give performances close to that which 
was anticipated. These checks performed on the machine give, 
in addition to rather precise dynamic characteristics for the 
coefficients selected, a good indication of the effect of certain 
other variations in selected coefficients. This technique will be 


illustrated in Example 5. 
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CHAPTER 5 


MODIFICATIONS INSTALLED ON THE FOURIER SYNTHESIZER 


The considerable labor required in cross plotting from the 
x and y tapes to the P (z) plane may be eliminated by the use ofa 
complex - vector adder connected to x and y signal leads. A 
servo driven two-axis recorder was usedforthis purpose. Timed 
pulses were introduced on the x and/or y axis to provide an in- 
dication of selectedvalues of 6. The harmonic synthesizer sweeps 
out from 0° to 360° at a uniform rate and @ is therefore propor- 


tional to time. 


A more satisfactory method of observing and recording the 
polynomial contour (P (z) plane) plot was found by utilizing a 
standard cathode ray oscilloscope in place of the two-axis recorder, 
and a Polaroid Camera for recording. The associated circuitry 
and additional components required to modify the harmonic syn- 
thesizer to provide a suitable computer are described in Appen- 
dix A. This resulted in a much smoother and more accurate plot 
since the mechanical and servo response difficulties inherent in 
the two axis recorder were eliminated. Also, since the area of 
primary interest on the P (z) plane is in the vicinity of the x axis 
cross-over point (where y =0), this area can be enlarged to cover 
the entire scope. Great latitude is possible in scaling both the 
synthesizer and scope, permitting substantial flexibility in select- 
ing the regions to be more closely examined. Through the use of 
the Polariod camera and multi-exposure, as many runs as neces- 
sary can be superimposed on one record for comparison. The 
sweep speed of the trace is slow enough so that it is practicable to 
follow the spot with a grease pencil directly on the face of the 
scope or on tracing paper. Thus the camera may not be required 
unless a permanent accurate record is required. 


o3 


This superposition of record technique is the only prac- 
tical method by which curves for constant @ and varying R can 
be obtained, since the synthesizer is inherently a constant R 


and variable @ device. 


The set of constant 0 curves are obtained merely by con - 
necting points of equal @ on the several constant R curves. For 
these records the synthesizer is not scaled for maximum amplitude 
for each run of different radii, but scaled as shown in Example 
o to provide the actual relationship of each radius plot to all the 


others. 


6 indication was provided by fabrication and installation 
of a cam on the synthesizer drive shaft as described in Appendix 
A. A microswitch actuated by the cam follower was installed on 
the synthesizer and the signal produced was used to blank out the 
scope for about AS , thus providing a @ indication. Cam teeth 
were cut at 8 = 90°, 120°, 135°, and 150°, and the cam was cut 
to blank out the scope from 180° to 360°, since the trace is 
symmetrical about the x axis and the last 180° would contribute 


nothing further. 


The x and y axes on the scope trace may be obtained simply 
by turning off the y and x scope amplifiers respectively. The 
x and y axes, being linear scales, may be calibrated in a manner 
identical to the method of calibrating the component tapes or as 
shown in Example 5. Thus, the x axis location of the constant 
term, (encirclement point) may be determined. The dynamic 
‘characteristics or roots of an equation with constant coefficients 
may be found to a satisfactory degree of accuracy quickly and 
easily without using large amounts of Sanborn Recorder tape. 
Various radii (R) runs are scaled and set into the synthesizer. 
Observing each trace for encirclements and the x axis cross - 
overs relative to the constant term on the x axis serves as a good 
estimate for the radius of the next run. Scaling and setting values 


into the machine requires very little time or effort. All roots of 
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the equation may be found inthis manner. A root is indicated 
when the trace passes through (or near, depending on the degree 
of accuracy desired) the encirclement point. The radius (R) 
used gives the W of the root, andthe 6 marks indicate the 
polar angle of the root, measured counter~clockwise from 7 

on the root plane | (ADR) = 71 - 6| , as well as the damping ratio 


[¢ = cos (ADR)| . No Sanborn. tape need be used unless 
desired. When a more precise value of (ADR) is required than 
that obtainable from the five @ indications provided by the cam, 
a Single tape can be run for the root value of the radius R and 


a more accurate value of @ obtained. 


For higher order equations, fifth or higher, it may be 
desirable to enlarge the scope presentation, particularly when 
searching for roots at the smaller radii, so that the immediate 
vicinity of the encirclement point occupies the entire presentation. 
When the origin no longer appears at the far right side of the scope, 
as with an enlarged presentation, it may be necessary to re - run 
previously scaled plots where the x cross overs are known in 
order to re~scale the portion of the negative real axis that ap- 
pears on the presentation. With some higher order equations it 
may be advantageous to divide through by the constant term prior 
to making any runs, thus shifting the origin to the constant term 
on the negative real axis. Thus the origin becomes the encircle- 


ment point and calibrating the x axis is no longer a problem. 


The principles, procedures and equipment just discussed 
are particularly well adaptable to equations with fixed coefficients 
and an adjustable constant term. Problems of this type arise 
frequently, particularly in control system and servo problems 
involving an adjustable open loop sensitivity. Runs are made as 
previously discussed for a few radii and the 0, and (ADR) at 
x cross-overs are plotted on the root plane. The value of the 
x axis intercept is the constant term for the equations thus formed 


with roots at the WO and ¢€ indicated by the cross-over. Or- 
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dinarily just a few such runs, with their corresponding roots and 
value of the constant term will be sufficient to fair in a reason- 
ably accurate root locus on the root plane plot. The effect on the 
magnitude and (ADR) of the roots of varying the constant term 
(or open loop gain) may be determined from the tape or photo 


record. 


The analysis of the case where all the coefficients and 
the constant term are adjustable will be illustrated by an example 
in the next chapter. This case arises when system parameter 
Specifications are to be set which will enable the resultant system 


to meet the required dynamic characteristics. 
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CHAPTER 6 


ILLUSTRATION OF THE ANALYSIS TECHNIQUE 


EXAMPLE 5: Adjustable Coefficient Case 


A typical performance equation of a lead modified positional 


servomechanism is 


1 


2 
Pp + Cy P 


where I 


+ 


— 
— 


l+yTp = Bigey Fp i 
1+7TP S(ps)| Aom l+7T p 5(ps) Ain M intf) (1) 


equivalent moment of inertia, referred to the 


controlled member 


damping coefficient, referred to the controlled 


member 


performance function of the lead modifier net~ 


work, where y <l. 


characteristic time of the lead modifier 


sensitivity of the positional servomechanism 


for angle in, torque out. 


interference torque 
angular position of the controlled member 


command input angle 


Laplace Transform operator 


Oo” 














Cc DE 
If we define: cae = 7 = 4 2 —(Ps) 
oy vy S 
ep) e P ane. == — . 
n 
"(ps ) (ps) 
Equation (l) becomes 
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The characteristic equation is: 
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ee es +(e +1) 'p+ 


l 
— = 0 4 
PARE aeons 
ee j9 Le : 
let p=z=Re ? 20 ips) a, nr Dee Cc 
b_ 3 2 
Then P (z) - roe a (a+b) z + (ab+l) z (5) 


or P(R,6)~2 = R® e398? + (atb) R® e 1°94 (abtl) Re I?) 6) 
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Taking reasonably broad limitations based on servo principles as 


a starting point for the analysis: 


.58 Sips) < 15 O65 7 <5 Bis oy 


Hence: l<a< 3 Web < 3 weeec< | 
and ,38 ab <9 


Therefore: 1.3<a+b< 6 Les ab sal 16 . IS ema 


Values of the second coefficient, a+b, were picked which 


covered the range of variation of atb; i.e. a+tb=1.3, 3.5, and 6. 


However, choosing a value for a+b restricts the possible 
variation of the third coefficient, ab+l, since b must remaina 


real number physically. 


Thus, if atb = 1. 3, 


a=1,3-b 
and abtl=x = (1.3-b) b+l 
Bebo = liseoeet 


Transposing and solving for b, 


Ueales VAS) (7) 
2 


bes 


The maximum value for x for which b remains a real number 
is obtained by setting the quantity under the radical equal to zero 


and solving for x: 


4 (x-1) = 1.69 
_ 1.69 nee 
a a +1 = 1.42 
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Therefore, lo 3 <abrl= 142 when atb=l. 3 


The value of b is determined from Equation (7) by the 


substitution x = 1. 42. 


Similarly, the permissible variations of ab+l when atb=3.5 


and a+b=6 were determined and found to be 
1.3 < abtl< 4. 06 when atb = 3.5 


and 1.3 < ab+1< 10 when atb = 6 


On the basis of the preceding estimates, runs were made 


for various R's with the following adjustable coefficient settings: 


Run —atb ab+l 2 

ea 1.3 1.3 ashes 

1-21 1.3 1. 42 655 2< 6.5 
2-1 3.5 1.3 3 <2<.9 
2-2 3.5 4.06 L.75< 2s 17.5 
3-1 6 1.3 3528.5 
3-2 6 10 gS 2 Ot 


These runs were photographed and are presented on the 
succeeding pages. The conclusions drawn from the examination 
of the photographs appear with the records. 
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Figures 13, 14, 15, 16, 17 and 18 For Example 5 Follow. 
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For any general case where it is desired to select coeffi - 
cients Specifying component parameters, the effect on the dynamic 
characteristics of varying any coefficient may be readily observed. 
This is shown in the following series of photographs, Figs. 19 
and 20. 


The equation used is of the form: 


P(z) = BORINAE i iS a: = 0 


A was taken as 1.3, 3.5 and 6. 

B was taken as 1.3, 10 andl19Q. 

C is unrestricted and to be selected 
from whatever specifications are 


imposed. 


It will be observed that this is simply an extension of the 
estimated boundaries of coefficients of the previous specific 
problem. Some of the same runs were used to demonstrate that 
there is nothing unique about any one set of runs, They might 
be equally applicable to any one of a variety of problems. 


For the purposes of illustration, one set of runs will be 
analyzed for possible conclusions relating the dynamic charac- 
teristics and the coefficients. All other sets for any problem 


are handled ina Similar manner. 


Photograph 3-3 representing 
Eze 2? 4 62° +19z + C = O 


is sketched in Fig. 21 for clarity although constant @ curves may 


be drawn directly on the photograph if desired. 
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Figures 19 and 20 which follow show the effect of the 


variation of the equation coefficients. 
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X= -400 





Analysis from ‘ 
Constant R curves a 


10<C 215 lcomplex pair SSW <6 
lrealroot .52a<1.0 


15<C <20 1 complex pair Aw SS 
1 real root l<a<3 


252C<65 1 complex pair 3<w, <4 
1 real root 3<a<5 


652C 2140 1 complex pair Ae =o 
1 real root 65 <a<4 


i “140 110 -80 -65 -40 - 





P(z)=z5+Az*+Bz+C 


10-5 -2 





0-20 45 


~\ 
‘\ 


‘\ 
\ 


a @=120° 
\ ADR=90° 
Neo 
\ C 
\ 
\ 





0<C<20 1 


— 


C= 20 


— ps 


20<C <30 1 
1 


30<C 240 1 
1 


2.5 
“Ty oe ; 


Analysis from 
Constant @ curves 


ve” 135° 
ADR=45° 





va ADR=60° 


complex pair 30° <ADR<45° 
real root ADR= 0 
complex pair ADR = 30°w = 4) 
real root (a # 1,5) 


complex pair 3P< ADR <45° 
real root 


complex pair 45°<ADR<60° 


real root 


402C 2250 1 complex pair 60°<ADR<90° 


250 <C 


1 real root 


1 complex pair ADR > 90° 
1 real root 


Unstable since one complex pair 


occurs in right half plane 


Fig. 21 Polynomial Contours from Photographs 3-3 Constant R and Constant 6 Curves. 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS: 


In this thesis, the mapping theorem was used in estimating 
the dynamic characteristics of physical systems from the char- 
acteristic equation of the performance function. Modification of 
a Fourier Harmonic Synthesizer greatly facilitated the application 
of the theorem, and enabledthe polynomial contour to be obtained 
directly. The presentation of the polynomial contour on a stand- 
ard cathode ray oscilloscope made it possible to apply the encir- 


clement criteria conveniently and with considerable accuracy. 


It was found that input signal levels too low to be indicated 
by the Sanborn Recorder pens still produced accurate polynomial 


contours on the scope presentation. 


The difference between the computed polynomial contour 
check points and the corresponding points on the photo records 


was not measurable. 


By expanding the scope presentation, any portion of tne 
polynomial contour may be closely examined. This procedure is 
useful when investigating critical, or limiting values of the charac- 


teristic equation coefficients. 


Examples have been included which indicate that the tech- 
niques presented are applicable either to polynomials with fixed 
coefficients or to polynomials whose coefficients are adjustable. 
In the case of adjustable coefficients, the method presented pro- 
vides for the possibility of setting the limits on the coefficients 
which will allow the roots of the system characteristic equation 
to remain within a specified sector on the root complex plane. 


13 


This will allow the dynamic characteristics of the system to meet 
the prescribed specifications. The performance functions of the 
system components are also determined and are represented by 
by the coefficients. This knowledge would be very helpful to the 


preliminary designer. 


Although the examples presented here are primarily chosen 
from servo problems, the technique is equally applicable to air- 
craft stability problems or to other problems where the dynamics 


determined by the characteristic equation are of interest. 


RECOMMENDATIONS: 


For the purpose of illustrating the analysis techniques for 
the adjustable coefficient case, a typical problem involving an 
equation of the third degree was chosen (Example 5). It is recom- 
mended that these procedures be applied to actual physical prob- 
lems of fourth degree or higher. The results of this investigation 
indicate that the techniques developed apply equally well to equa- 
tions of higher order without excessively increasing the amount 


of computation required. 


It is further recommended that the feasibility of inverting 
this technique be studied. The shape of the polynomial contour 
is readily adjustable by the synthesizer potentiometers while the 
machine is in operation. Once the relation between the polynomial 
plane contour and the root plane is fully appreciated, it should be 
possible to adjust the synthesizer to produce polynomial contours 
representing the desired dynamic characteristics. By starting 
with the potentiometer settings, and inverting the scaling tech- 
nique, the equation coefficients providing the specified dynamic 


characteristics may be determined. 


If such a technique can be developed, it should facilitate 
further the problem of estimation of the dynamic characteristics 


of physical systems. 
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APPENDIX A 


IMPLEMENTATION OF THE HARMONIC SYNTHESIZER 
MODIFICATIONS 


The harmonic synthesizer used in obtaining data for this 
thesis was constructed by the Instrumentation Laboratory of the 
Massachusetts Institute of Technology. The complete circuit 
diagrams of the synthesizer are available at the Laboratory’) The 
mechanics of the synthesizer operation will be outlined in order to 


clarify the reasons for the modifications made. 


This synthesizer is composed of the fundamental and 23 
harmonic components. Figure A = 1 illustrates the operation of 
the fundamental component of the synthesizer. An electric motor, 
through a gear train, furnishes the drive for the rotation of the 
rotor coils, which are 90° apart electrically. The moving coils 
are excited by a 400 ~ voltage in the stator coils. The signal vol- 
tages are picked off at the summing resisters, amplified, and sent 
to a Sanborn two-channel recorder. Since the rotors for each 
component are 90° apart electrically, the resulting signal voltages 
will have a phase difference of 90° and thus correspond to cos 6 
and sin 6, respectively, as the motor turns the rotor shaft through 
one complete revolution of 360°. Multiplication by R is done by 
the potentiometers. Thus, the output voltages of the fundamental 


component in the synthesizer correspond to Ry cos @ and R, sin 6, 


l 
or to the real and imaginary terms, respectively, of 


e Ve Ry cos @+] Ry sin @ 


nd 


Similarly, if the rotor coils of the2 ~ harmonic component 


have twice the angular velocity of the fundamental, the output 
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voltages of this component correspond to 


aes R, ec 379 = R, cos 20+ j Ry Sitiezee 
1) 36 , 
anal = Re J = Ro cos 36+] R. sin 30, etc. 
Bt R, = BR, Ro ~ ee and R3 = ee then the sum of the 


first three components of the synthesizer is 


2 


P (R, 6) = (R° cos 30+ AR“ cos 20+BR os 6) 


+j (R°sin 3 6+ AR 2 sin26+ BR sin 6). 


Thus, the real and imaginary parts of any complex harmonic 
function (within the limits of the machine) may be readily obtained 
by expressing it as a polynominal in polar form, such as 


6 n-l , j(n-1) 6 


a Sa 
P(R,6)-a,=R"e +a,_,R +. 


, ta, Re j2e 


+a, Re 19 
The characteristic equation of any physical system may be 


expressed in this manner. 


The real and imaginary parts of the function~- henceforth 
referred to as the x and y signals, respectively - were applied 
to the x and y inputs of a standard Dumont cathode ray oscilloscope 


which was fitted with a Polaroid Land Camera 


A z- dot generator was designed in order to provide scope 


illumination at the moments when the 400 ~ excitation signal was 


16 


at its maximum, as shown in Figure A - 2, 


These pulses 
i appear on the face 
of the scope at the 
} | 7 rate of 400/ second 
\ | and appear to be a 
| ) a steady dot, which 
: traces the complex- 


vector sum of the 

x and y signal vol- 
Figure A-2: Signal Voltage tages. Thus, a 
with Z4-Dot Generation visual represen- 

tation of the poly- 

nominal contour is 

presented, which is 

the first requirement for the application of the encirclement cri - 


terion. 


Finally, a cam-actuated microswitch was installed on the 
rotor shaft of the fundamental component, which provided timing 
indications corresponding to certain values of 6. The final 180° 
of the dot travel, being symmetrical to the first 180°, was blanked 
out. Figure A -3 shows the Fourier Harmonic Synthesizer after 
modification, and Figures A-4,5, and 6 are the circuit diagrams 
of the x and y signal amplifiers and the z-dot generator. Figure 7 


is a sketch of the cam. 


In operation, as the dot traced out the plot of the character- 
istic equation polynomial , the microswitch cam follower opened 
contacts, thus blanking the scope, duringthe indicated intervals. 
The scope camera, set for time exposure, photographed the trace 


fora permanent record. 
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APPENDIX B 


Tabulation of Numerical Valucs Used in the Determination 
of the Effect of Varying the Coefficients in the Characteristic 


Fquation of Example 4: 


P(z) = 2° + (atb) 22 + 2+b 
=o +lat2t wi) 2° +(2b wi ato) a4 a geet 
where atb = ea + 20. Ww. 
yo 1 i o'gs cea: z 
- GO, @ ve 
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APPENDIX C 


Tabulation of Fourier Harmonic Synthesizer Scale Settings Used 
in Example 5. 


The basic equation used in this example is 
P(z) = 2° + (atb) 2° + (abtl) z+2 
or, with z= Re ive 
P(R,0)-2 = R°e 394 (atb)R* e 19 + (abt) Re J 
NG, ol a eM ele je 
el 2 14+40 


1 : 
=e = lle al 


and the real part of P (R, 8) 4 = Ee = ae + (a+b) Re + (abil )ake 
@=0 


Similarly, at 0=7, e JF e a Ji = -14+j0 


zg leg e 14". aan Le) 
ele qe jst Tee 
and the real part of P (R, @) ~ : ae -R°+(atb) R~ (ab+l)R 
0=7 


Since, in this example, eee = 6, the largest potentiometer 
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setting, S3 will correspond. ta R° 


mx, 


a throughout. That is, 


3 
max 


3 


division of all numerical coefficients by R = 6 = 216 gives 


the scale setting for the potentiometer in each run. 


These tabulations are carried out in the following tables, 


as indicated. 


Sample calculation: 


R=6 R= 36 R® = 2i6 

1.3R7 = 1,3 (36) = 46.8 

13R = 13(6) = 7.8 
a = (216) = ous F = 216+ 46.8 + 7.8 = 270.6 
3° 216 216 . 0 3 3 3 
So = ste (1-3R°) = == (46.8)=0.217 F =-216 + 46.8 - 7.8=-177 
2 - aie‘ DIG ee rt . ! 


S 


| = 51g (1. 3R) = ag (7.8) = 0.0361 
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